A kinetic equation for linear fractional stable 
motion with applications to space plasma 

physics. 



Nicholas W. Watkins and Dan Credgington ' 

Physical Sciences Division, British Antarctic Survey (NERC), Cambridge, 

CBS OET, UK 

Raul Sanchez 

Fusion Energy Division, Oak Ridge National Laboratory, Oak Ridge 

TN 37830, USA. 

Sandra C. Chapman 

Centre for Fusion Space and Astrophysics, University of Warwick, Coventry 

CV4 7AL, UK 



Abstract 

Levy flights and fractional Brownian motion (fBm) have become exemplars of the 
heavy tailed jumps and long-ranged memory seen in space physics and elsewhere. 
Natural time series frequently combine both effects, and Linear Fractional Stable 
Motion (LFSM) is a model process of this type, combining alpha-stable jumps with 
a memory kernel. In contrast complex physical spatiotemporal diffusion processes 
where both the above effects compete-dubbed "ambivalent" by Brockmann et al. 
(2006)-have for many years been modelled using the fully fractional (FF) kinetic 
equation for the continuous time random walk (CTRW), with power laws in the 
pdfs of both jump size and waiting time. We derive the analogous kinetic equation 
for LFSM and show that it has a diffusion coefficient with a power law in time 
rather than having a fractional time derivative like the CTRW. We develop earlier 
comments by Lutz (2001) on how fBm differs from its fractional time process coun- 
terpart. We go on to argue more physically why LFSM and the FFCTRW might 
indeed be expected to differ, and discuss some preliminary results on the scaling of 
burst sizes" and "durations" in LFSM time series, with applications to modelling 
existing observations in space physics. 

Key words: Plasma transport, anomalous diffusion, time series modelling 



Preprint submitted to Elsevier Science 



20 March 2008 



1 Introduction: The need for non-Brownian models of complexity 
in space plasma physics. 



Fractional equations have seen at least two applications in space physics. 
They have been proposed in fra ctional kinetics-based mode ls of microphys- 
ical processes in space plasmas (jMilovanov &: Zelenyil . l200ll ). because of the 
link between fractional ki netics and observed anomalous transport in magnet- 



ically confined plasmas (jBalescul . Il995l ; ISanchez et all l2005t iBalescu et al. 



20051 ) They are also used to model magnetospheric activity, because they can 
describe the non-Gaussian probability density functions and long-range tem- 
poral memory seen in real time series, su ch as tho s e from auroral indices and 
some in-situ solar wind quantities (e.g. IWatkins It is not obvious, 

however, that it must be the same type of fractional equation in both these 
contexts-i.e. reversible microphysical plasma transport equations versus equa- 
tions for the evolving features of time series from a macroscopic variable. In 
this paper we investigate th is open question , prom pted in part by the stimu- 
lating suggestions made by IZaslavsky et al.l (120071 ) in response to our earlier 



paper (IWatkins et al.l . 120051 ) 



A historical preamble (section 2) will introduce the two main non-Brownian 
properties seen both in spatiotemporal anomalous diffusion and in many 1- 
dimensional time series. They are heavy tailed pdfs and long-range memory, 
Mandelbrot's "Noah" and "Joseph" effects respectively. We will also introduce 
some of the models in use to study them, notably the Levy flight (ordinary 
Levy motion or oLm), fractional Brownian motion (fBm) and the continuous 
time random walk (CTRW). We will then discuss the highly topical problem 
of "ambivalent" processes where the Noah and Joseph effects compete. 

Having described the phenomena and the main types of model, in section 
3 we will recap the main diffusion-like equations modifying Wiener Brownian 
motion (WBm) which have embodied these models in the parallel literatures 
on stochastic time series modelling and anomalous diffusion. While we are 
of course aware that concepts and models flow back and forth between these 
overlapping communities we feel using such a structure has made it clear how 
the models have developed, and exposed some surprising lacunae. 

We note that the same fractional diffusion equation holds for t he pu r e Noa h 
effect in both the oLm and CTRW descriptions. As noted by iLutzl (120011 ). 
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such an identity no longer holds in the case of the corresponding paradigmatic 
models for the pure Joseph effect. We show these two different equations- for 
the fractional time process and fBm-and discuss why they differ. Prompted by 
this difference we go on to look at two corresponding self-similar models which 
unite the Noah and Joseph descriptions, one described by the fully fractional 
CTRW and one by linear fractional stable motion (LFSM). This codification 
of the set of equations allows us to show that, surprisingly, a kinetic equation 
is "missing" from the literature, that for LFSM. We give a simple derivation 
for it by direct differentiation using its well-known characteristic function. 



In section 4 we explore the application of LFSM to modelling the "burst" 



measu red on magnetospheric and solar wind time series (e.g. 
( 2000h N ). We make simple scaling arguments building on result of 



Freeman et al. 




Kearnev & Maiumdar 



ent SOC" behaviour (c.f. IWatkind (120021 )) and compare them with numerics. 



In section 5 we discuss the interesting fact that LFSM shows additive rather 
than rational ambivalent behaviour, i.e. the self-similarity exponent H is an 
additive rather than rational expression. 



In section 6 we summarise and then conclu de by considering the i mplica- 
tions of this work for the arguments made by IZaslavsky et al.l (120071 ) which 
advocate the use of the fully fractional CTRW in space physics time series 
modelling. 



2 Time series and anomalous diffusion: Phenomena and models 



2.1 Non-Brownian time series 



2.1.1 Self similarity, the Hurst effect and H. 



Hurst's observation (the "Hurst effect") of the anomalous rate of growth of 
range in hydrological time series, such as the height of the river Nile, was one of 
the first natural phenomena for which the need f or a non-Brownian description 
was recognised. iMandelbrot fc Van Nesd (119681 ) explained the Hurst effect as 
being due to long range dependence in time, which they referred to as the 
"Joseph effect". They encapsulated the Joseph effect in their seminal model, 
fractional Brownian motion (fBm). Like WBm, fBm has the property of self 
similarity under a dilation in time where At goes to XAt: 



x(XAt) = \ H x(At) 
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Throughout this paper we follow lEmbrechts fc Maejimal (120021 ) in denoting 
by H the self similarity exponent defined this way. Unlike WBm where H is 
always 1/2, in fBm H takes values between and 1. 



2.1.2 Long range dependence, R/s and the Joseph exponent J. 



To describe the growth of rescaled range (R/ s) due to persistence. iMandelbrot fc Van Ness 
( 119681 ) used the Joseph exponent J, where 



R/s ~ r 



(2) 



Normal diffusion and time series modelled by random walks have J = 1/2 = 
H. For the particular case of fBm we still have H = J, and indeed the absolute 
value of the displacement x and the mean square of displacement < x 2 > grow 
with time as J and 2 J respectively. Because of this equivalence, in the case of 
fBm, one can use Mandelbrot's "R/s" method to measure H. More generally, 



however, as he emphasises (e.g lMandelbrotl . 120021 . p. 157) H and J will differ 



and R/s measures the latter. We hence do not use the potentially confusing 
term "Hurst exponent" in this papei0- 



2.1.3 Heavy tails, the Noah effect and the stability exponent fi. 

A second type of non-Brownian phenomenon had also been recognised by 
Mandelbrotl (119631 ). This was the non-Gaussian increments, with heavy power- 
law tails, 



P[x) 



x 



-(i+aO 



(3) 



seen in financial time series and also in many natural ones. In contrast to the 
Joseph effect he called this the "Noah effect" . He proposed a second paradig- 
matic model, ordinary Levy motion (oLm), for cases when the anomalous 
behaviour of the time series originates entirely from this effect, rather than 
long temporal memory. In oLm the index \i runs from to 2. 

A self- similarity exponent H remains defined for oLm, by H = l/fi. How- 
ever it was recogn i sed early on that J would present conceptual subtleties 
(jMandelbrotl . |2002| ; IMandelbrot fc Wallisl . 119691 ) stemming from the fact that 
in oLm all moments of order greater than \x are infinite, including the second 
order moment used in the R/s method. A finite data series drawn from an 



2 In the particular case of fBm we also have that (3 = 2H + 1 
is the exponent from the power spectral density (S(f) ~ f~^) 



2 J + 1, where -j3 
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infinite variance process thus presents "pseudo- Gaussian" behaviour to diag- 
nostics such as R/s or the variance. In particular J takes the value it would 
do for an equivalent process with Gaussian amplitude distribution. J is 2 for 
both oLm and ordin a ry Br ownian motion, for example. As pointed out by 



Mandelbrot fc Wallid (119691 ) this is an advantage if one actually just wants 



to measure the degree of temporal memory and not any other source of self- 
similarity, but it also means that R/S can longer be used to measure th e full 
self-similarity exponent H (see also the discussion in Mandelbrot ( 20021 )). 



2.1-4 Ambivalent time series: Noah and Joseph in competition 



Real time se ries usually do not ex h ibit j ust one or the other of these two 
limiting cases. Mandelbrot fc Wallisl (119691 ) thus proposed that the effects 
modelled by fBm and oLm could be combined in a more general self simi- 
lar additive model, "fractional hyper bolic" motion, now refe r red to as linear 
fractional stab le motion, LFSM, (e.g. lEmbrechts &: Maejimal . 120021 ) . Mandel- 
brot (p. Ill in IShlesinger et all (119951 )) described LFSM later as "one of the 
"bridges ... combining the infinite variance feature ... and the global depen- 
dence feature He had evidently found LFSM unsatisfactory as a financial 
model, and remarked in 1995 that it "[had not] found concrete ... use". It 
has ho wever, by now, been applied to p roblems as diverse as comm unications 
traffic (ILaskin et all 120021). geophysic s (IPainter fc Patterson! I1994T). magneto - 
spheric physics (IWatkins et all 120051 ) and solar flares (IBurnecki et all 120081 ) . 



In this selfsimilar and stable random walk there are two contributors to H, 
i.e. H = H(L, J) = L + J ' — 1/2. This is reminiscent of the physical finding that 
in nonequilibrium statistical mechanics the spatial and temporal correlation 
lengths need longer be the same 7^ £11). As we might expect the self- 
similarity exponent of the Noah effect L depends on the spatial exponent fi 
controlling the amplitude of the steps in the walk, via L = l//x. The temporal 
exponent (3 controls the Joseph exponent J for long range memory by J = 
(3/2 — 1/2. In the specific case of finite variance processes such as fBm we can 
equate (3 to the power spectral exponent. 



2.2 Modelling anomalous diffusion 



All the above time series models have counterparts describing the now 
widel y-recognised natura l phenomenon of anomalous (non-Brownian) diffu- 
sion (IKlafter et allll996l ) in space and time. Particularly relevant to the space 
plasma case is the physics of non-equ i librium sys tems, and tur b ulent diffu- 
sion, see e.g. the reviews of IZaslavskyi (120041 ) and iBalescu et al.l (120051 ). The 
paradigmatic model here has been the continuous time random walk (CTRW), 
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where the object studied is the joint probability P(x, t) that a random walker 
makes a jump of size x after waiting for a time t. The f actor ising form where 
P(x,t) = P(x)P(t) has been particularly well explored, with power laws in 
P(x) and P(t ) being used to model lar ge jumps and long range memory effects, 
respectively (iMetzler fc Klafterl . lioooh . 



Diffusive systems in which the Noah and Joseph effects compete have re - 
cently become highly topical in complexity science ( Brockmann et all 2006 ). 
These authors used the CTRW with power laws in both P(x) and P(t) to 
model the space-time dynamics of dollar bills carried by travelle rs. They coined 



the ph rase "ambivalent" diffusion for such a "mixed" process. IWatkins et al. 



( 120051 ) had earlier drawn attention to the need for this type of model for 
space physics time seri e s app lications, but proposed the use of LFSM in such 



cases. 



Zaslavsky et al.l (120071 ) criticised their approach. While agreeing that 



an ambivalent model was indeed needed IZaslavsky et al.l (120071 ) argued that it 
should be of the CTRW type and that a procedure was thus needed to define 
Levy jumps in such natural t ime series. Very recent ly LFSM has been used to 



model solar flare time series (IBurnecki et al.l . 120081 ). 



3 Non-Brownian equations of motion codified 



3. 1 Random walks, the Central Limit Theorem and the diffusion equation 



The physical phenomenon of Brownian motion and its mathematical ide- 
alisation as the Wiener process are central to equilibrium statistical physics. 
Their links to the Central Limit Theorem (CLT), and the problem of limit 
distributions in general, are thus key elements in the relation of mathemat- 
ics to physics. Furthermore their embodiment in diffusion equations has been 
essential both to practical applications and to physical understanding. 

The familiar form of the diffusion equation is 

d d 2 

-k;Pwbm(x, t) = D——P WBM (x, t) (4) 



Parallel developments in the mathematical modelling of time series, in 
nonequilibrium diffusion, and turbulence, have all allowed progress in describ- 
ing systems for which a description beyond Brownian motion is needed. We will 
here focus mainly on the first two, and on additive models where increments 
are added to produce a random walk. Multiplicative models are particularly 
natural in turbulence but our discussion of these will be limited to the issue 
of multifractality in section 5. 
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3.2 Heavy tailed jumps (the Noah effect) 



There are many possible ways of modifying the CLT. One involves retaining 
self similarity, but allowing long-range correlations in space or time. Such ex- 
tensions have been studied in at least two formalisms. One is the mathematical 
theor y of stable processes ( Embrechts fc Maejima . 20021 ; Samorodnitskv and Taqqu . 
1994) and other i s the continuous time rando m walk (CTRW) ( IMetzler fc Klafterl . 



2000 



Zaslavskyl . l2004j ; iBalescu et all 120051 ). Both formalisms arise when the 



finite variance assumption of CLT is relaxed, and both lea d to a description by 
a diffusion equation with a fra ctional derivative in space (IPaul &: Baschnagell . 
19991 ; IMetzler k Klafterl . I2004J ) 



d 

■^PoLM(x,t) 



D—P 

OX 11 



OLM 



X,t) 



(5) 



The two formalisms must be equivalent here because they both correspond 
to iid infinite variance processes (Levy flights) and so must be equivalent 
manifestations of the extended CLT. 



3.3 Long ranged temporal memory (the Joseph effect) 



The situation is more subtle when the iid assumption is relaxed, rather than 
the finite variance one. Relaxing independence is one way to break iid. One 
wa y to relax independen c e is i ntroducing; long range dependence, as studied 
by iMandelbrot fc Wallisl (119691 ) . They used a self-affme process with a mem- 
ory kernel, and named it fractional Brownian motion. In the CTRW formal- 
i sm it was instead introdu ced via a power law distribution of waiting times 
( Lindenberg fc West . 19861) s o was known as the fractional time (or temporal) 
process (FTP, see also iLuta (120011 )). 



3.3.1 Long ranged memory in the CTRW picture: the Fractional Time Pro 
cess 

The modificat i on to the diffusion equation (4) that corresp onds to the FTP 
■alal 



3DC 

is ( iBalakrishnanl . Il986l : IMetzler fc Klafterl . 120001 : iLutd . 120011 ) 



d 2 

dt ^P FT p(x,t) = D—P FTP (x,t) 



(6) 



Note that we do not include the term describing the power law decay of 
the initial value here or in subsequent CTRW equations (it is retained and 
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discussed by Metzler fc Klafter ( 2000l ). see their equation 40). The fractional 
derivative in time corresponds physically to a power law in waiting time be- 
tween jumps. The notation f3 just indicates that this exponent need not nec- 
essarily be related to the (3 we use in the following. In all the following cases D 
is no longer the Brownian diffusion constant but simply ensures dimensional 
correctness in a given equation. 



3.3.2 Long ranged memory from a self-similar process: fractional Brownian 
motion 



Contradictory statements exist in the literature concerning the equivalent 
kinetic equation for fBm corresponding to equation (6) for FTP. It has some- 
times been asserted ( IZaslavskyl . 12004 IWatkins et al.U2005l ) that (6) is also the 
equation of fBm. However th e solution Pftp of (6) is now known to be of a 
non-Gaussian form, given by (IMetzler fc Klafterl. 120001) in terms o f Fox func- 



tions ; whereas the pdf of fBm is by definition (IMandelbrotl . Il982l ; iMcCauley 



2004 ) Gaussian but with the standard deviation "stretching" with time as t H . 
The correct diffusion equation for fBm must thus, as noted by Lutz (2001), 
be local in time: 



dt 



FBM 



2Ht 2H ' 1 D 



dx 2 



FBM 



(7) 



given, to our knowledge, firstPl by Wang fc Lung ( 1990 ). It can be seen by 
trial solution to have a solution of the required form . Furth er properties of the 
solutio ns to equation (7) have been given by iLuta (120011 ) . In addition [Lutz 
( 120011 ) has given a corresponding fractional Langevin equation. fBm and the 
FTP are very different in the sense that, although both include temporal corre- 
lations, the kinetic equation for the former is Markovian in time and the latter 
is not. This is in spite of both processes being non-Marko vian in respe ctively 
the Wiener and CTRW senses! This underlines (see also Lutz ( 200ll )) that 
these two closely related ways of modifying the CLT by introducing temporal 
dependence doin fact have strikingly different structures. 



With hindsight we may understand why equations (7) and (6) differ on 
physical grounds. The fBm approach (7) seems to us to be more macroscopic 
in spirit, in that a time dependent diffusion coefficient is being imposed versus 
the more microscopic approach of prescribing a pdf for waiting time. How- 
ever we would an t icipat e that the nonlinear shot noise formalism studied by 
Eliazar fc Klafterl (120061 ) might allow one to derive one formalism as a limit 
of the other. 



note the diffusion coefficient was given in Feder ( 19881 ) 
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3.4 Heavy-tailed jumps and long range memory together: two possible diffusion- 
like equations for ambivalent processes 



3-4-1 Fully fractional C TRW and Brockmann et al's (rational) ambivalence: 

Similar questions have been asked over the years about the natural gen- 
eralisation of equation (6) for the fractional time process to allow for Levy 
distributions of jump lengths as well as power-law distributed waiting times. 
The resulting equation (corresponding in particular to Brockmann et al's "am- 
bivalent process" ) is fractional in both in space and time: 

-q^-Pffctrw{x^) — D-^P FFC TRw{%it) (8) 



Again, as with the FTP the solution for this process is known (IKolokoftsov et al. 



20011 ) not to be a stable (or stretched stable) distribution but rather a convo- 



lution of such distributions. 

3-4-2 Additive ambivalence: LFSM 

As with ffim, equation (7), there is a process, linear fractional stable motion 
(LFSM), which generalises the one represented by eqn (7) to the ambivalent 
case. 

Its pdf Plfsm can be defined via its characteristic function 
Plfsm= I e ikx eM-a\k\^ H ) (9) 



see for example ( Laskin et al. . 20021 )). 



We see that LFSM has a Levy-like characteristic function: exp(— a^^t^ ) 
where the effect of /i no longer being equal to 1/H is for the effective width 
parameter to grow like t^ H . The characteristic function has the correct fBm 
limit, where we take fi as 2. For fBm at any given t the characteristic function 
is a Gaussian with width growing as t 2H = t^ -1 = t 2J . 

That LFSM is a general stable self-affine process can be seen by taking 
k' = kr H whereby we obtain 



LFSM 



r H Mx/t H ) (io) 



This is a stable distribution of index fi and a prefactor which ensures the 
self-similarity in time discussed in section 2.1.1 
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The kinetic equation satisfied by (9) which generalises (7) can be found by 
direct differentiation of (9) with respect to time to give 



d 



g- t P LFSM = LiHat^- 1 J e ik *\krexp(-a\k\^ H ) (11) 



which, absorbing the constant a, and factors of 2n into D can be recognised 
as 

f) flu 

Q- t P LFSM = flHt^D — PLpsM (12) 



using one of the standard definitions ( Paul fc Baschnagel . 1999 . p. 110) of a 



fractional derivative d^/dx^. Surprisingly equation (12) seems not have been 
given before in either the physics or mathematics literature. 

The appropriate limits may be easily checked; in particular \x = 2 gives the 
equation of fBm found by Wang and Lung. We also remark t hat LFSM should 



be a sp ecial case of the nonlinear shot noise process studied by lEliazar fc Klafter 



(120061 ) which may allow further generalisation of the equation we have pre- 



sented. 



4 Scaling properties of ambivalent processes and their relatives 



4-1 Burst durations in LFSM 



Watkind (120021 ) reviews the use of power law pdfs in "burst" size and du- 



ration derived from time series to infer the presence of self- organised crit- 
icality (SOC) in the magnetosphere and solar wind (e.g. by iFreeman et al. 



(120001 )). Qualitatively such behaviour may also simply be a n artifact of a 



self sim ilar (or multifractal) time series (a possibi lity raised bylFreeman et al. 
( 2000 )). To clarify this we have elsewhere (e.g. Watkins ( 20021 )) advocated 
the testing of SOC diagnostics using controllable self similar models. In this 
section we present an example (Figure 1) of the pdf of burst duration for a 
simulation of LFSM. LFSM was simulated using a direct FFT-based method 
( Chechkin &: Gonchar . 2000l ; Watkins et all 2005 ). Burst duration was defined 
in the standard manner as being the period for which a time series exceeds a 
given threshold. \x and f3 were chosen as 1.5152 and 1.58, giving a subdiffusive 
H v alue of 0.45. These parameters are qu ite typical of the geomagnetic index 
AL dWatkinsl . 120021 ; IWatkins et aT1 . l2005h . 
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pdf D(T) of burst durations T for fl_m fi=1 .51 52 ,p=1 .58 



* Simulation 
Predicted from 2- H- 1.55 



x ^ v 

* N s. 

* V . 

* 



io' 2 - 



10"' 10° 10 1 10 2 10 3 10 4 10 5 10 s 10 ? 10' 



Fig. 1. Probability density D(T) of a burst of duration T for simulated LFSM 

The expected behaviour of burst duration for LFSM is not completely ob- 
vious from the literature. We have overlaid a power law of exponent 2 — H 
as this gives the scaling of the isoset (set of threshold crossings) of a general 
fractal and should, we believe, thus be independent of the detailed nature of 
the model. This preliminary comparison suggests reasonable agreement. We 
should note however that detailed agreement with measured exponents is not 
attempted, and we would not necessarily expect it as LFSM is a very over- 
simplified model. 
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4-2 Burst sizes in LFSM 



Figure 2 shows a representative pdf D(e) of burst size e, again for simulated 
LFSM. Burst size was defined as being the area above the threshold while the 
time series exceeds it. fi and j3 were chosen as before. The power law scaling 
is less straightforward and there appears to be a turnover for small e. 



We now develop a simple scaling argument, following iKearney fc Majumdar 
( 120051 ) to predict the scaling of the tail of the pdf of LFSM for larg e e. It is 
thus o ne candidate toy model fo r such "apparent SOC" behaviour (c.f. IWatkins 
( 120021 )). IKearney fc Majumdarl (120051 ) considered the zero-drift Wiener Brow- 
nian motion (WBm) case. Rather than their full analytic treatment we first 
recap their heuristic argument for a burst size (area) A using first-passage 
time tf. This may then be adapted to isosets and then to LFSM. They first 
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note that for WBm y(t) ~ t 1 / 2 for large t. Then, defining A by 



A = J y(t')dt' 



(13) 



the integration implies that large A scales as tj . Simple inversion of this 
expression implies that tf must scale as tf ~ A 2 / 3 . We independently have the 
standard result for first passage time for WBm: P(t f ) ~ t~ 3/2 . To get P(t f ) 
as a function of A i.e. P(tf(A)) one needs to insert the expression for tf as a 
function of A in ab ove equation, and in addition will need a Jacobian. After 
these manipulations Kearney fc Majumdarl (120051 ) found that 



P(A) ~ A" 4/3 



(14) 



In the zero-drift but non-Brownian case we can argue that y(t) ~ t H for 
large t, but rather than first passage times we define an isoset-based burst size 
e using the set of up and down crossing times {ti} such that the integral from 
U to t i+1 is 

e = J y(t')dt' (15) 



As remarked by iFreeman et al.l (120001 ) the points {ti} form an isoset of the 



fractal curve. For a fractal curve of self similarity exponent H and dimension 
D = 2 — H they have dimension 1 — H. In consequence the probability of 
crossings over a time interval r goes as t 1 ~ h giving an inter-event probability 
scaling like t~^~ h \ The pdf for inter-event intervals in the isoset thus scales as 
r -{2-H) ^ fas same scaling as for the first passage distribution in the Brownian 
case. 

The rest of the argument goes as before. If y ~ t H then e ~ t\ +H so tj ~ 
e -( 1 +- ff ). Folding in the scaling of the inter-isoset intervals (ij), P(ti) ~ tj H ^ 
we can get P(ti) as a function of e (i.e. P(ti(e))) by inserting an expression 
for £j as a function of e in the above equation, and again using a Jacobian. 
We find: 

P(e) ~ e' 2 /( 1+H ) (16) 



which we can check by taking the Brownian case where H = 1/2 to retrieve 
P(e) ~ e -4 / 3 . Comparing the numerical trial in Figure 2 with this predicted 
scaling exponent of —2/(1 + if) is again encouraging, but detailed comparison 
with data is postponed to a later paper. 
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pdf D(e) of burst sizes e for fLm u=1 .51 52 ,p=1 .58 



x Simulation 

— Predicted from 2/(1+H)=1 3793 



to" 2 * : x » » 



10 10 10 10 10 10 10 10 10 10 

e 

Fig. 2. Probability density D{e) of a burst of size e for simulated LFSM 
5 Two sorts of ambivalence: additive and rational 

5.1 Ambivalent processes 



Despite the primacy of mathematical descriptions, well-chosen verbal nomen 
clature has always been very important to physics both in guiding thought 
and in condensing prev ious insights. A recent example is the coinage by 
Brockmann et al.l (120061 ) of the phrase "ambivalent process" . We believe that 
the term should in fact be extended to all processes in which long ranged tem- 
poral correlation and long-ranged amplitude jumps compet e, not just the fac- 



torised (fully fractional CTRW) diffusive process studied by lBrockmann et al. 



(120061 ). Their coinage prompted us to recognise that LFSM is also an am- 
bivalent process and we feel that this wider usage would encourage further 
research on the inter-relationships between these two types of ambivalent pro- 
cess and others. We feel it also has the merit of focusing further attention on 
how LFSM and the CTRW differ from Levy Walks and multifractals. 



5.2 Rational versus additive ambivalence 



Examining equation (12) it is clear that LFSM is also an ambivalent process 
in Brockmann et al's sense, but equally clear that it differ s subtly from the 
CTRW . The type of ambivalence previously identified by iBrockmann et al. 
(120061 ). comes from the factorised FF CTRW. We may rename it rational be- 
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cause of clearly different scaling exponents for the average spatial displacement 
< x(t) >~ t^ 1 - 1 compared to LFSM. Rather than £(1) being given by a ratio 
(3' / fi of temporal and spatial exponents, as happens for the FF CTRW, for 
LFSM it is an additive function of (H = J + L — 1/2) of the relevant analogous 
exponents or functions thereof: 

H = [l/p] + [(3/2 - 1/2] - 1/2 (17) 



The difference in scaling seems to arise because the LFSM path itself is defined 
as a convolution (when the long time limit is taken) between a memory kernel 
and a stochastic jump, and thus factorises in Fourier space. Conversely, for 
the decoupled CTRW, it is the jump pdf which factorises into a temporal and 
spatial term. 



5.3 A coupled space-time model: the Levy walk 



As opposed to such models which are often in some sense factorised or con- 
volved by hand, the underlying physics itself may lead to space and time being 
explicitly coupled. A good example is the Levy walk which was introduced to 
deal with the problem of infinite variance by directly coupling the distribution 
of flight times to a long-tailed jump distribution. Well-known a pplications to 



natura l systems have included the Swinney group's experiments (IKlafter et al 



( 119961 )) on particle tracer transport in turbulence. In a Levy walk the exponent 
for < X 2 (t) > has several possible functional dependencies on the control pa- 
rameters, a clear diff erence from the CTRW a lready remarked on by Metzler 



and Klafter (p. 30 of iMetzler fc Klafterl (120001 )). It is of course no coincidence 



that early evidence for Levy walks came from turbulence. Space and time 
are unavoidably coupled in multiplicative processes, and Levy walks describe 
diffusion on such fields. 



5.4 Beyond self- similarity: the multifractal 



Other more complicated fractal models such as multifractals have also been 
extensively applied. Although only one of these four broad classes of model 
may be applicable to a specific case, availability of a wider choice of scaling 
exponents should enable a better discrimination of the underlying mechanism 
in at least some complex systems. 
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6 Conclusions 



In this paper we studied the question of whether one would expect the same 
equation to describe a time series as an anomalous diffusive pr ocess, prompted 
both by the sti mulating; suggestions of IZaslavsky et al. ( 2007 ) and the highly 
topical work of iBrockmann et al. I (l2006h on "ambivalent' diffusion processes 
where the Noah and Joseph effects compete. A codification of diffusion-like 
equations showed that a kinetic equation was "missing" from the literature; 
the one corresponding to LFSM. We gave a simple derivation for it by direct 
differentiation of the well-known characteristic function of LFSM. We then 
made a preliminary exploration of how LFSM could model the "burst" sizes 
and durati ons previously measu red on magnetospheric and solar wind time 
series (e.g. iFreeman et al.l (120001)). We made sim ple scaling arguments build- 
ing on a result of Kearney fc Majumdarl (120051 ) to show how LFSM could 
be one candidate explanation for such "apparent SOC" behaviour and made 
preliminary comparison with numerics. We also discussed the interesting fact 
that LFSM shows additive rather than rational ambivalent behaviour, i.e. the 
self-similarity exponent H is an additive rather than rational expression. 



All this has consequences for the interesting arguments of IZaslavsky et al 
(120071 ). We emphasise that, perhaps surprisingly, one cannot not simply equate 
the limiting cases (i.e. FTP and FF CTRW) of the CTRW with their fflm 
and L FSM analogues. This point had earlier been made for fBm by iLutz 
(120011 ); we have shown here that the problems he pinpointed must also apply to 
LFSM. They have interesting consequences for the applicability of the CTRW 
to persistent time series-derived data, leading us to believe that models which 
were designed with time series applications in mind, such as LFSM may be 
more suitable in many cases. 



Future work will include testing the above conclusions with different non- 
Fourier based generators for LFSM. We also plan to conside r other stochas- 



tic p rocesses, both FARIMA and nonlinear shot noises (c.f. (IBurnecki et al. 



20081 )). to allow generalisation of the above initial investigations into burst size 
and duration. The prevalence of natural processes showing heavy tails and/or 
long ranged persistence suggests a relevance well beyond the initial area of 
application in space physics. 
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